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ABSTRACT 


We present the measured Sunyaev-Zel’dovich (SZ) flux from 474 optically-selected MaxBCG clusters that 
fall within the Atacama Cosmology Telescope (ACT) Equatorial survey region. The ACT Equatorial region 
used in this analysis covers 510 square degrees and overlaps Stripe 82 of the Sloan Digital Sky Survey. We 
also present the measured SZ flux stacked on 52 X-ray-sclected MCXC clusters that fall within the ACT 
Equatorial region and an ACT Southern survey region covering 455 square degrees. We find that the measured 
SZ flux from the X-ray-selected clusters is consistent with expectations. However, we find that the measured 
SZ flux from the optically-selected clusters is both significantly lower than expectations and lower than the 
recovered SZ flux measured by the Planck satellite. Since we find a lower recovered SZ signal than Planck , 
we investigate the possibility that there is a significant offset between the optically-selected brightest cluster 
galaxies (BCGs) and the SZ centers, to which ACT is more sensitive due to its finer resolution. Such offsets 
can arise due to either an intrinsic physical separation between the BCG and the center of the gas concentration 
or from misidentification of the cluster BCG. We find that the entire discrepancy for both ACT and Planck can 
be explained by assuming that the BCGs are offset from the SZ maxima with a uniform random distribution 
between 0 and 1 .5 Mpc. In contrast, the physical separation between BCGs and X-ray peaks for an X-ray- 
selected subsample of MaxBCG clusters shows a much narrower distribution that peaks within 0.2 Mpc. We 
conclude that while offsets between BCGs and SZ peaks may be an important component in explaining the 
discrepancy, it is likely that a combination of factors is responsible for the ACT and Planck measurements. 
Several effects that can lower the SZ signal equally for both ACT and Planck , but not explain the difference in 
measured signals, include a larger percentage of false detections in the MaxBCG sample, a lower normalization 
of the mass-richness relation, radio or infrared galaxy contamination of the SZ flux, and a low intrinsic SZ 
signal. In the latter two cases, the effects would need to be preferentially more significant in the optically- 
selected MaxBCG sample than in the MCXC X-ray sample. 


Subject headings: cosmic microwave background - galaxies: clusters: general - galaxies: clusters: intracluster 
medium 
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1. INTRODUCTION 

Galaxy cluster properties may follow simple scaling laws 
reflecting their self-similarity. This possibility has given 
credence to their use as cosmological probes. Cosmolog- 
ical parameters have been obtained by X-ray, optical, and 
most recently SZ surveys of clusters (e.g., Vikhlinin et al. 
2009; Mantz et al. 2010; Rozo et al. 2010; Sehgal et al. 201 1; 
Benson etal. 2011). At the same time the scaling laws that 
feed ir.to these parameter constraints continue to undergo 
scrutiny. 

Recent millimeter-wavelength data have opened a new 
window whereby these scaling relations can be robustly 
checked against SZ flux measurements. The SZ cluster sig- 
nal has been predicted to have a low-scatter correlation with 
cluster mass (e.g., Moti et al. 2005; Nagai 2006). If true, 
this would make SZ-detected clusters an excellent tracer of 
structure growth in the Universe (e.g., Wang & Stcinhardt 
1998; Haiman et al. 2001; Holder et al. 2001; Caristrom et al. 
2002). As steps towards understanding the SZ-mass rela- 
tion, several studies have shown that the SZ cluster sig- 
nal correlates well with X-ray signals (e.g., Bonamentc et al. 
2008, 2011; Andersson et al. 2011), dynamically determined 
masses (Sifonetal. 2012), and weak-lensing determined 
masses (e.g., Marrone et al. 2011). 

In particular, the Planck satellite recently reported a 
good agreement between the measured and expected SZ- 
mass relation for a sample of X-ray-selected clusters 
(Planck Collaboration et al. 2011a,b). However, a similar 
comparison for optically-selected clusters yielded an am- 
plitude of SZ flux lower than expected by about a fac- 
tor of two, with an even larger discrepancy for lower-mass 
clusters (Planck Collaboration et al. 2011c). An analysis by 
Draper et al. (2011) using data from the WMAP satellite found 
a similar result, however with larger uncertainty. Hand et al. 
(201 1). using data from the Atacama Cosmology Telescope 
(ACT) (Swetz et al. 2011) and stacking luminous red galax- 
ies, also suggested a low SZ flux for optically-selected halos. 
Among the possible explanations could be that either the SZ 
signal is not a robust tracer of galaxy clusters and groups or 
that optical selection techniques are somehow biased. 

Here, we investigate this discrepancy by stacking optically- 
selected clusters in millimeter- wavelength data from ACT that 
overlaps Stripe 82 of the Sloan Digital Sky Survey (York et al. 
2000). We also measure the SZ flux for X-ray-selected clus- 
ters as a consistency check. Understanding these scaling rela- 
tions will have important implications for cluster astrophysics 
as well as for their use in cosmological studies. 

This paper is organized as follows. Section 2 discusses the 
data sets used in this analysis. Section 3 describes the method 
used to measure cluster SZ flux. Results are presented in Sec- 
tions 4 through 6 and discussed in Section 7. 

2. DATA SETS 

Below we describe the catalogs of optically-selected and X- 
ray-selected clusters and the millimeter-wavelength data used 
to measure the cluster SZ fluxes. We note that throughout this 
work Msooc refers to the mass within R 5000 which is the ra- 
dius within which the average density equals 500 times the 
critical density of the Universe at the cluster redshift. Sim- 
ilarly, Mzoom gives the mass within Rzoom, the radius within 
which the average density equals 200 times the mean mat- 
ter density of the Universe at the cluster redshift. A fidu- 
cial cosmology of = 0.27, ft a = 0.73, and h =■ 0.71 is 


also adopted (Komatsu et al. 2011), where H(z) = HqE(z) = 
(h x 100 km s" 1 Mpc ~ l )E[z) and E(z) = [O m (l +z) 3 + 0 A ] I/2 . 

2.1. The MaxBCG Optical Cluster Catalog 

The MaxBCG Optical Cluster Catalog consists of 13,823 
clusters selected from Data Release 5 (DR5) of the Sloan Dig- 
ital Sky Survey (Koester et al. 2007a,b). The clusters were 
selected from a 7500 deg 2 area of sky using the observation 
that cluster galaxies tend to be the brightest galaxies at a given 
redshift, share a similar red color, and are spatially clustered. 
The catalog consists of clusters that fall in the redshift range 
of 0.1 < z < 0.3 and have a richness measure, A^oo/n* within 
10 < Nzoom < 190. The richness is defined as the number of 
red- sequence galaxies with L > 0 AL * (in the i band) within 
a projected radius of /? 200 m. The catalog provides the BCG 
position (RA and DEC), photometric redshift, richness, BCG 
luminosity, and total luminosity of each cluster. Applying the 
cluster detection method to mock catalogs suggests that the 
catalog should be 90% pure and 85% complete. 

Mass estimates of the clusters in the MaxBCG sample 
were derived by Sheldon et al. (2009) and Mandelbaum et al. 
(2008a) using weak gravitational lensing. Johnston et al. 
(2007a) and Rozo etal. (2009) used those mass determi- 
nations to construct richness-mass (A^oo™ - Mooc) relations. 
Rozo etal. (2009), in particular, used the masses derived 
from Mandelbaum et al. (2008a) due in part to the au- 
thors’ careful treatment of photometric redshift uncertainties 
(Mandelbaum et al. 2008b). Rozo et al. (2009) also stacked 
the MaxBCG cluster catalog on X-ray maps from the ROSAT 
All-Sky Survey (Voges et al. 1999) and used the Lx -M re- 
lation from Vikhlinin et al. (2009) as a prior to inform their 
richness-mass relation. Thus the richness-mass relation of 
Rozo et al. (2009) is expected to be consistent with an L x -M 
relation from X-ray clusters. Rozo et al. (2010) applied this 
richness-mass relation to the MaxBCG sample of optically- 
selected clusters and found a cosmological constraint on erg of 
0“8(£WO.25) 041 = 0.832 ± 0.033 assuming a flat ACDM cos- 
mology. Throughout this work we define the A^oom “Mooc re- 
lation as given by Eqs. 4, A20, and A21 of Rozo et al. (2009). 

2.2. The MCXC X-ray Cluster Catalog 

The Meta-Catalog of X-ray detected Clusters of galaxies 
(MCXC) is presented in Piffaretti et al. (201 1). The MCXC 
cluster catalog is based on publicly available data from a num- 
ber of different X-ray catalogs including the ROSAT All-Sky 
Survey and comprises 1743 clusters. The catalog provides 
the position (RA and DEC), redshift, X-ray 0.1 -2.4 keV 
band luminosity (Tsooc), mass (Mooc), and radius (J? 5 ooc) of 
each system. The redshift distribution of this catalog goes 
from about 0.05 to 1. The L x - M relation derived from the 
MCXC clusters in Piffaretti et al. (2011) is consistent with 
that of Pratt et al. (2009) and Vikhlinin et al. (2009). The 
Vikhlinin et al. (2009) L x -M relation was used to derive a 
cosmology constraint on cr 8 from a sample of X-ray clusters 
that is a subsample of the MCXC catalog. From this analy- 
sis they found cr 8 (O m / 0.25) 047 = 0.813 ±0.013 (stat) ±0.024 
(sys) (Vikhlinin et al. 2009). Other authors have found sim- 
ilar constraints on 0 % from ROSAT and other X-ray cluster 
samples (Henry et al. 2009; Mantz et al. 2010). 

2.3. Millimeter-wave Data from the Planck Satellite 

Given that the optically-selected MaxBCG cluster catalog 
and the X-ray-selected MCXC cluster catalog yield consistent 
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Lx -Mrelations and constraints on one would expect both 
cluster samples to yield consistent Tsooc-^sooc relations. 22 
However, in a set of papers presented by the Planck collabo- 
ration (Planck Collaboration et al. 201 lb,c) it was found that 
the Ysoq c -Msqqc relation for X-ray-selected clusters from the 
MCXC sample agrees with expectations, whereas the normal- 
ization of the 7500 c ~ M 500 c relation for the optically-selected 
MaxBCG sample was lower than expectations by about a fac- 
tor of two. For both cases, expectations were based on X-ray 
derived cluster profiles from Amaud et al. (2010). For the lat- 
ter case, the expectation also folded in the A^oo/n - Mooc rela- 
tion of Rozo et al. (2009). The Planck data used in the above 
analysis consists of six HFI channel millimeter-wave tem- 
perature maps as described in Planck HFI Core Team et al. 
(2011). This data set comprises the first ten months of the 
survey and covers the full sky. 

2.4. Millimeter-wave Data from the Atacama Cosmology 
Telescope 

The Atacama Cosmology Telescope is a six-meter tele- 
scope operating at an altitude of 5200 meters in the Atacama 
Desert of Chile. The telescope site allows ACT to observe in 
both the northern and southern hemispheres. In this work, we 
use millimeter-wave maps covering two regions of sky: one 
spanning 510 deg 2 over the celestial equator and one span- 
ning 455 deg 2 in the southern hemisphere. The Equatorial 
region consists of a 4.5° -wide strip centered at a declination 
of 0° and running from 20*20'” through 0* to 03*50”. The 
Southern region consists of a 7° wide strip centered on -53° 
and extending from 00*12” to 7*10”. Both sky regions were 
observed over the 2008, 2009, and 2010 observing seasons 
at 148 and 218 GHz. The Equatorial region overlaps the 
Sloan Digital Sky Survey (SDSS) Stripe 82 and thus over- 
laps 492 clusters in the MaxBCG catalog. The Equatorial 
plus Southern regions combined overlap 74 clusters in the 
MCXC catalog. For a more detailed description of the ACT 
instrument, observations, and data reduction see Fowler et al. 
(2007); Swetz et al. (2011); Marriage et al. (2011); Das et al. 
(2011); Hajian et al. (2011); Dunner et al. (2012). 

3. MEASUREMENTS OF SZ FLUX 

3.1. Multi-frequency Matched Filter 

We use a multi-frequency matched filter to extract 
the thermal SZ signal from clusters as described in 
Haehnelt & Tegmark (1996) and Melin et al. (2006). The fil- 
ter in Fourier space is given by 

*00- [P(k)]- 1 • r(k), (1) 

where r(k) has the components 

r i ,(k) = r / (k) 7 ,^(k). (2) 

Here j v is the frequency dependence of the thermal SZ sig- 
nal for frequency u, r'(k) is the profile of the cluster in 
Fourier space, and B v { k) is the profile of the instrument beam 
in Fourier space. P(k) is the power spectrum of the noise, 
both astrophysical and instrumental. The astrophysical noise 
sources for cluster detection include the primary lensed mi- 
crowave background, radio galaxies, infrared galaxies, Galac- 

22 1500c is the SZ flux within /?5oo c> dividing out the frequency dependence 
of the SZ signal. 


tic emission, and the SZ background from unresolved clus- 
ters, groups, and the intergalactic medium. Since the power 
from the cluster thermal SZ signal is subdominant to these 
astrophysical sources (as evidenced by Luekeret al. (2010); 
Hall et al. (2010); Fowler et al. (2010); Das et al. (201 1)), we 
approximate the power spectrum of the total noise as the 
power spectrum of the data itself. Here 

A = Jd 2 k [r(k)]' ■ [P(k)]- 1 • [r(k)]J _1 (3) 

is the normalization of the filter that ensures an unbiased esti- 
mate of the cluster signal. 

3.2. SZ Model Template 

We use for the filter’s spatial template the empirical 
universal pressure profile of Amaud etal. (2010) derived 
from X-ray observations of the REXCESS cluster sample 
(Bohringer et al. 2007). The three-dimensional pressure pro- 
file is given by 

where x = r/r s , r s = Rsooc/csoo, C 500 = 1 .156, a. = 1 .0620, /3 = 
5.4807, and 7 = 0.3292. The normalization of this profile is 
arbitrary for the purposes of the matched filter. We essentially 
measure this normalization for each cluster when we apply 
this filter to our maps. The SZ signal is given by the projected 
gas pressure, so we describe the filter template by integrating 
the three-dimensional profile above along the line of sight. 
Thus 

P 2D {d) = f “ 2 P iD ( ^P+8 2 D a (z) 2 ) dl , (5) 

Jo 

where l max = 5R$qo c andZ^(z) is the angular diameter distance. 
The filter is truncated at 5R 5 oo c /D A (z) = 50 5O oc , which contains 
over 95% of the signal, as was done in Melin et al. (201 1); 
Planck Collaboration et al. (2011b,c). 

\ 

3.3. Application of Filter 

Before filtering the maps, we establish uniform noise prop- 
erties by creating an effective weight map that has pixel-wise 
effective weights given by w eff = (“ + ^) _1 where w t is the 
pixel weight for the fth frequency. Here weight is defined as 
the number of observations per pixel normalized by the obser- 
vations per pixel in the deepest part of the map. We multiply 
the 1 48 and 2 1 8 GHz ACT data maps pixel-wise by the square 
root of the effective weight map. After we apply the filter to 
create a filtered map, we divide the filtered map pixel-wise by 
the square root of the effective weight map. 

We apply the matched filter above following the procedure 
given in Planck Collaboration et al. (2011c). For each cluster 
in the MaxBCG catalog that falls in die ACT coverage region, 
we create a unique matched filter using the Njoom - M 5 oq c re- 
lation of Rozo et al. (2009) to determine each cluster’s Msooc 
and subsequently R 5 m c from the N 2 oo m value given in the cat- 
alog. We also derive D A (z) for each cluster from its photo- 
metric redshift. When the filter is applied to the map, the 
pixel coincident with the location of file cluster center in the 
filtered map should have a value equal to the normalization 
of the two-dimensional SZ template given by Eq. 5. To sim- 
plify extraction of the desired quantity, we normalize P 10 {0) 
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FlG. 1 . — Recovered F 500 values from 447 simulated clusters embedded in 
ACT 148 and 218 GHz maps (black circles). These clusters were simulated 
to match the properties of the clusters in the MaxBCG catalog (Koester et al. 
2007a) that overlap the ACT equatorial region. The simulated clusters were 
placed at random locations within the ACT maps. The input F 500 values arc 
shown as blue squares. 
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FIG. 3. — Black circles show recovered X 500 values from 446 embedded 
simulated clusters where the assumed cluster centers used for 7 50 o recovery 
are offset from the input centers with a random distribution given by the black 
solid line in Figure 2. "Mcent" stands for miscentered. The input values are 
shown as blue squares. 


itself to equal unity when integrated over 0 in two dimensions 
from zero to 5 # 500 ^ Thus the pixel value recovered at the 
cluster center position after applying the filter is 7^^ . Here 

= y srL,> where WL is ^ integrated projected 
SZ signal within a cylinder of radius 5Rsooc We use a geomet- 
ric factor of = (0.986/1 .8 14))^^ given in Appendix A 



Fig. 2. — The black solid line shows the distribution of offsets between 
BCG and X-ray gas peak from 208 clusters found in both the MaxBCG and 
MCXC cluster catalogs. The red dashed line shows the distribution for the 
subsample of rich clusters (N 200 m > 35). The blue dotted line shows the 
same for the subsample of poor clusters (A^oom < 35). There are 100 and 108 
clusters in each subsample respectively. 


of Melin et al. (2011) to convert from 7 5 ^ qoc to Y^ c , the in- 
tegrated SZ flux within a sphere of radius Rsqqc- Through- 
out this work we plot f 50 o = Y^qE~ 2 ^(z)(D a (z)/ 500 Mpc) 2 , 
where 7 50 o = 7^* , as in Planck Collaboration et al. (201 lc). 

4. SIMULATED ACT SZ SIGNALS 
4.1. Embedding Simulated Clusters in ACT Maps 

In order to test the analysis pipeline which applies the SZ 
extraction procedure discussed above, we use simulated clus- 
ters embedded within the ACT data maps at random locations. 
Using the information in the MaxBCG catalog for the 492 
clusters that fall within the ACT Equatorial region, we create 
a unique SZ profile for each cluster using Eq. 5 above and the 
cluster Rsooc and z given in the catalog. We then add each 
simulated SZ cluster to the 148 and 21 8 GHz ACT Equatorial 
maps, placing it at a random location and scaling the thermal 
SZ signal to give it the appropriate frequency dependence in 
each map. The simulated SZ signal is also convolved with the 
appropriate ACT beam prior to embedding it within each ACT 
map. We then exclude any simulated clusters from further 
analysis that happen to be within 5' of a point source detected 
at greater than 5a in either the 148 or 218 GHz maps. We also 
exclude any clusters that are within 10' from the edge of the 
map. As an additional cut, we exclude all clusters that are in 
noisy parts of the map where the local value of the effective 
weight map is less than 15% of the maximum value. These 
cuts leave 447 clusters. Similarly below, when we extract the 
SZ signal from the real cluster positions using the MaxBCG 
and MCXC catalogs, we apply the same cuts discussed above. 

In Figure 1, we show the results of this SZ extraction pro- 
cedure. Here the simulated clusters are binned by the opti- 
cal richness given in the MaxBCG catalog. The blue squares 
show the input model, and the black circles show the recov- 
ered signal of the simulated clusters embedded at random lo- 
cations within the ACT maps. The error bars are given by 
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FIG. 4. — Measured J 500 values for 52 MCXC X-ray-selected clusters 
(PifEaretti et al. 201 1) that fell within the ACT equatorial and southern sur- 
vey regions (black circles). Also shown are expected ^500 values based on 
measured cluster X-ray properties (blue squares). A cluster profile model 
from Arr aud et al. (2010) was assumed for determining both measured and 
expected F 500 values. 


Fig. 5. — Measured K 500 values for 474 MaxBCG optically-selected clus- 
ters that fell within the ACT equatorial survey region (black circles). Ex- 
pected T 500 values are shown as blue squares. Both measured and ex- 
pected values assume the N 20 Qm - Afsooc relation from Rozo et al. (2009) 
and the Amaud et al. (2010) cluster profile. Red triangles are the measured 
values from the Planck satellite for a sample of 13,104 MaxBCG clusters 
(Planck Collaboration et al. 2011c). 


a/y/N where a 2 is the variance of 7soo in each richness bin, 
and N u the number of clusters in each richness bin. The vari- 
ance dominates the uncertainty in each bin, as demonstrated 
in Figure 4 of Planck Collaboration et al. (201 lc). 23 There is 
good agreement between input and recovered signals. 

4.2. Effect of Cluster Miscentering 

The exercise above shows we should expect excellent re- 
covery of the SZ flux for clusters given in the MaxBCG cata- 
log provided the N 2 00 m - M 500 c and M 5 qq c - Y S qo c relations are 
correct, and the cluster properties (position, redshift, # 200 ™) 
listed in the catalog arc accurate. However, one source of un- 
certainty identified in Johnston etal. (2007b) (section 4.3) is 
the positional accuracy of the cluster center. In the MaxBCG 
catalog, the cluster position is given by the location of the 
Brightest Cluster Galaxy (BCG). Johnston et al. (2007b) sug- 
gested rwo reasons why the BCG found by the MaxBCG clus- 
ter identification algorithm may be offset from the true clus- 
ter center. One is that the true BCG may be intrinsically 
offset from the dark matter center or center of the gas con- 
centration presumably due to unrelaxed behavior (e.g., merg- 
ers). Another reason may be that the BCG could be misiden- 
tified by the cluster finder. Johnston et al. (2007b) explore 
this latter effect with mock optical cluster catalogs and find 
that a richness-dependent fraction of the clusters have accu- 
rately identified BCGs while the rest are miscentered due to 
BCG misidentification following a Rayleigh distribution with 
a scale parameter cr equal to <7™ sc = 0.42/j _, Mpc. The dis- 
tribution of the intrinsic BCG offset from the gas center is 
unknown. 24 

23 Since no scatter in the A^oo« - Msooc relation has been included here, 
these errer bars reflect only the SZ flux recovery error. 

24 Note that this offset docs not arise from pointing uncertainties in optical. 
X-ray, or millimeter-wave instruments. It is due to either BCG misidentifica- 
tion or cluster astrophysics. 


We explore the potential offset between BCG and gas cen- 
ter by studying the subset of clusters in common to both the 
MaxBCG and MCXC catalogs. From the full catalogs (which 
have 13,823 and 1743 clusters respectively), we identify 208 
clusters that are in common in both catalogs. Here we define 
a cluster as matched in both catalogs when the identified clus- 
ter redshifts are within A z < 0.015 and the projected physi- 
cal separation of the identified cluster centers is less than 1.5 
Mpc. Using these 208 clusters, we plot the fraction of clus- 
ters as a function of separation between BCG and X-ray peak 
in Figure 2. The solid black line shows the offset distribu- 
tion for the 208 clusters. Half of the clusters have offsets less 
than 0.1 Mpc, while the other half have a roughly flat offset 
distribution between 0.1 and 1.5 Mpc. 25 

The 208 clusters were then divided into rich clusters (with 
# 200 m > 35) and poor clusters (with N 2 oo m < 35). A richness 
cut of # 200 m = 35 divides the 208 clusters into roughly even 
subsamples of ~ 100 clusters each. The dashed red line in 
Figure 2 shows the offset distribution for the rich clusters, and 
the dotted blue line shows the distribution for the poor clus- 
ters. Similar distributions are found for both subsamples, with 
the poor clusters showing slightly more offset. 

Using the simulations described above, we explore the ef- 
fect of cluster miscentering on SZ flux recovery. We use the 
same simulated clusters embedded at random positions in the 
148 and 218 GHz ACT maps as before. However, when re- 
covering the SZ fluxes, we use positions for the clusters that 
differ from the true cluster gas centers with a random distri- 

25 We note that the choice of 1 .5 Mpc is somewhat arbitrary. When we 
allow matches within 1 Mpc, we find 189 clusters. We also find that if we 
allow matches within 3 Mpc, Figure 2 plateaus instead of dropping to zero 
at large separation. The pLateau is due to poor clusters, as the distribution 
of rich ones does tend to zero. The size of a cluster (/? 200 m) is usually less 
than 2 Mpc, so this suggests that some of these poor cluster matches may be 
spurious. 
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Fig. 6. — The black open circles show the measured Tsoo values from simu- 
lations with an offset between identified cluster center and gas concentration 
center given by a uniform random distribution between 0 and 1.5 Mpc. The 
solid black circles, blue squares, and red.triangles are the same as in Figure 
5. The values from simulations have been shifted in the x-axis for clarity. 

bution given by the black solid line in Figure 2, We also allow 
for 10% of the clusters to be false detections to match the pu- 
rity of the MaxBCG catalog (Koester etal. 2007b). Figure 3 
shows the result of the SZ recoveiy process. The black circles 
show the recovered Tsooc values for simulated clusters embed- 
ded in ACT data. To give a s,ense of the relation between Mpc 
and arcminutes, for the redshift range of the MaxBCG clus- 
ter sample (z e (0. 1 , 0.3)), 0.5 Mpc corresponds to about 2' to 
4.5'. Given the ACT beam size (1.4' at 148 GHz and 1' at 21 8 
GHz), we expect a decrease in the recovered 7>ooc. signal due 
to the amount of miscentering shown in Figure 3. How much 
the recovered SZ signal decreases depends on the noise prop- 
erties of the ACT data, which differ significantly from pure 
white noise. We discuss this further in Section 7. 

5. MEASURED ACT SZ SIGNALS 

5.1. Stacking X-ray Selected Clusters 

Within the 510 deg 2 ACT Equatorial region and the 455 
deg 2 ACT Southern region are located 74 clusters found in 
the MCXC catalog of X-ray-selected clusters (Piffaretti et al. 
2011). After making the cuts discussed in Section 4.1 to ex- 
clude clusters near bright point sources, near map edges, or in 
very noisy parts of the map, 52 MCXC clusters remain. Using 
the R 5000 M500 a and redshift of each cluster and the projected 
SZ profile given in Eq. 5, we calculate the expected mean 7 50 o c 
values in each Msooc bin shown as the blue squares in Figure 
4. The black circles in Figure 4 show the mean of the re- 
covered Tsooc values using the multi-frequency matched filter 
given in Eq. 1 and the projected SZ profile created uniquely 
for each cluster. The error bars are the error on the mean 
given by a /y/N, where o 2 is the variance and N is the num- 
ber of clusters in each bin. Figure 4 shows overall agreement 
between expected and recovered SZ signals for the X-ray- 
selected clusters. The reduced chi-squared is 0.76 using 7 de- 
grees of freedom. This is consistent with the agreement found 
by Planck Collaboration et al, (201 lb) for a larger sample of 


X-ray-selected clusters. 

5.2. Stacking Optically Selected Clusters 

In the ACT Equatorial region there are 492 MaxBCG clus- 
ters. This reduces to 474 clusters once the above men- 
tioned cuts are made. 26 Using the #200*2 - M S qq c rela- 
tion of Rozo et al. (2009) and the M 5 oo c - Ysooc relation of 
Amaud et al. (2010), we find the expected and recovered 7 50 oc 
values of the MaxBCG clusters using the method described in 
Section 3. Figure 5 shows the expected values as blue squares 
and the recovered values as black circles. The recovered sig- 
nal is significantly lower than the expected signal as well as 
the signal recovered by Planck from 1 3, 1 04 MaxBCG clusters 
(red triangles) (Planck Collaboration et al. 2011c). 

From Figure 3, we see that some amount of offset between 
BCG and gas peak can result in a lower measured signal than 
expected. However, this figure also shows that the offset dis- 
tribution given by Figure 2 does not result in a low enough 
measured signal to explain the measurement shown in Fig- 
ure 5. To investigate the amount of offset necessary to match 
the ACT measured 7soo values shown in Figure 5, we redo 
the analysis shown in Figure 3. However, this time we use 
an offset distribution that is uniformly random between 0 and 
1.5 Mpc. We also again allow for 10% false detections. The 
results are shown as the open black circles in Figure 6. The 
solid black circles, blue squares, and red triangles in Figure 6 
are the same as in Figure 5. We see that this amount of off- 
set between BCG and gas peak roughly matches the measured 
values. An extensive scan of offset distributions is beyond the 
scope of this work, but it may be that a more complex or re- 
fined distribution could give a better fit to the measurements. 

5.3. Stacking Optical Clusters Using BCG Dominant 
Subsample and New Richness Measure 

Using a subsample of clusters with " dominant BCGs" we 
examine whether the expected SZ signal is closer to the mea- 
sured values. Such a subsample may more closely correspond 
to an X-ray-selected subsample, and with such a subsample 
Planck Collaboration et al. (2011c) found better agreement 
between model and measurement. We follow the definition of 
"BCG dominant" used by Planck Collaboration et al. (2011c) 
which is defined relative to the quantity LscG/^tot-LBCG)* 
Here L^ x and Z, B cg are the R-band luminosities of the cluster 
and cluster BCG respectively. For a "BCG dominant" cluster, 
this ratio is larger than the average ratio for a given richness 
bin. From the sample of 474 MaxBCG clusters used above, 
126 are "BCG dominant". Figure 7a shows the recovered 7soo 
values versus the model expectation for this subsample. The 
Planck measurements of the sample of 1 3, 1 04 MaxBCG clus- 
ters (not a subsample) is included for reference. 

Recently a new measure of cluster richness was developed 
with less scatter than the measure presented in Koester et al. 
(2007a) (Rykoffetal. 2011). We test whether using this 
new richness measure will yield differing results with regard 
to measured versus expected SZ signal. A catalog with a 
new richness measure assigned to each MaxBCG cluster is 
available online. 27 While for the previous richness measure 
the cluster richness-mass relation was calibrated using weak 

26 Note that the ACT region is centered on Stripe 82, and extends beyond 
it. So fewer real clusters are cut than simulated clusters because the real 
clusters are only in Stripe 82, and thus are located towards the central part of 
the ACT map which has lower noise. 

27 http://kipac.stanford.edu/raaxbcg 
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Fig. 7. — Left (a): Shown are the recovered F 500 values as in Figure 5, except using a subsample of 126 "BCG dominant" clusters from the 474 MaxBCG 
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richness measure for the MaxBCG cluster sample given in Rykoff et al. (2011). Note that both the measured values (black circles) and expected values (blue 
squares) change using the new richness measure as compared to Figure 5. The red triangles are the Planck measured values using the old richness measure for 
reference. 


lensing, this new relation was calibrated with an abundance 
matching technique (Rykoff et al. 2011). Calibration of this 
new measure via weak lensing is still in progress. We use the 
richness-mass relation given in Eq. B6 of Rykoff et al. (201 1) 
to determine Msooc and subsequently R$oo c - Figure 7b shows 
the measured and expected I500 values using this new rich- 
ness measure for the 474 MaxBCG clusters. The red triangles 
show the Planck measurements using the old richness mea- 
sure for reference. Note that Planck would have different re- 
sults if they use the new richness measure. 

5 .4. Contamination from Infrared and Radio Galaxies 

To investigate whether infrared galaxies may be reducing 
the SZ decrement at the MaxBCG cluster positions, we re- 
cover the >500 values from the 474 MaxBCG clusters studied 
above using the single band 218 GHz ACT map alone. We 
compare that to >500 values extracted at 474 random positions 
within the 218 GHz map. 28 A positive correlation between 
MaxBCG clusters and infrared galaxies would result in neg- 
ative >500 values compared to the random sample. We find 
that both the MaxBCG cluster sample and the random sample 
have an SZ flux consistent with zero in the 2 1 8 GHz map, and 
we do not detect any significant excess of infrared signal cor- 
related with the MaxBCG sample. This is shown in Figure 8. 

We also cross-correlate the MaxBCG cluster catalog with 
the VLA FIRST catalog of radio sources to investigate how 
much radio galaxies may be reducing the measured SZ decre- 
ment. 29 The FIRST survey uses the NRAO Very Large Array 
and covers over 10,000 square degrees of sky to a sensitivity 
of about 1 mJy at 1 A GHz. This survey also overlaps with the 

28 For these measurements of Ksoo> we do not divide out the amplitude of 
the frequency dependence of the SZ signal, which is close to zero at the null 
frequency. So these are really measurements of -AT /r cm b. 

29 http://sundog.stsci.edu/first/catalogs/readme_l 2feb 1 6.html 


Sloan Digital Sky Survey. We cross-correlate to find the frac- 
tion of 474 MaxBCG clusters that have a radio source above 
a given flux threshold, within 5' of the identified cluster po- 
sition. Assuming a typical spectral index for radio sources of 
-0.7 (e.g., Condon 1984; Linetal. 2009), we choose a flux 
threshold cut of 50 mJy at 1 .4 GHz to yield sources above 
2 mJy at 1 50 GHz. Such a source convolved with the ACT 
beam would have a temperature increment of about 30 \iK at 
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Fig. 9 — Recovered F 500 values from 984 simulated clusters embedded 
in simuleted CMB maps at 148 and 218 GHz with Planck- like noise (green 
diamonds) or the corresponding level of white-noise (purple hexagons). The 
maps were convolved with ACT beams as well as a 5' Gaussian beam, so that 
they correspond to Planck resolution. The recovered positions are offset from 
the true positions with a uniform random distribution between 0 and 1 .5 Mpc. 
The input >500 values are shown as blue squares, and red triangles are the 
measured values from the Planck satellite for a sample of 13,104 MaxBCG 
clusters (Planck Collaboration et al. 201 lc). 


the white noise case in Figure 9 . The results of the Planck-likt 
case in Figure 9 are also very different from those of the sim- 
ulated ACT case shown by the open black circles in Figure 6. 
The two differences between the simulations are the beam and 
the noise. In the ACT data there is 1/f noise, atmospheric 
noise, and noise from the primary microwave background. 31 
All of this results in a redder noise spectrum than Planck's. 
The presence of red noise in the maps causes the matched 
filter to suppress more power on large scales than would be 
the case for white noise. This causes the filter to return a 
lower signal than actually exists if the signal is extracted from 
a position that is offset from the cluster center. If the signal 
is extracted at the cluster center, however, the matched filter 
will return the correct signal regardless of whether the noise 
is white or somewhat red (as can be inferred from Figure l). 32 

We demonstrate this effect in Figure 10, where the blue 
solid curve is a 5' Gaussian profile. This is roughly the fil- 
ter response for a cluster profile convolved with the Planck 
beam as Planck's noise is close to white (see Eq. 1). The red 
dashed curve shows a 1 .4' Gaussian profile, which had been 
Fourier transformed and had all Fourier modes set to zero for 
/ < 2000, roughly analogous to the effect of the ACT noise. 
The profile was then Fourier transformed back. One can see 
how a low flux recovery results from an offset from the center. 
The green dotted curve shows the same as the red curve except 
for a 5' Gaussian profile. The low signal from a miscentered 
position is still apparent. This is why merely smoothing the 
ACT maps to match the Planck beam would not allow us to 
reproduce the Planck measurement. 


148 GHz. Since the typical SZ signal from a cluster is about 
100 fiK to within a factor of a few when smoothed with the 
ACT beam, a 2 mJy radio source would start to significantly 
reduce the SZ decrement. We find that about 10% of the 474 
MaxBCG clusters investigated above have such a radio source 
within 5' of its identified center. This small correlation is not 
enough to explain the large discrepancy between measured 
and expected SZ signals shown in Figure 5, although it may 
be a contributing factor. 

6. SIMULATED PLANCK SZ SIGNALS 

To investigate what Planck's measured SZ signal would 
be if there existed the amount of offset modeled in Fig- 
ure 6, we again use simulations. We embed 492 simulated 
MaxBCG clusters in two sets of simulated CMB maps at 148 
and 218 GHz that have Planck-\ike instrument noise added 
to them. 30 We model the Planck noise using the the noise 
power spectra at 143 and 217 GHz shown in Figure 35 of 
Planck HFI Core Team et al. (2011). We allow for an offset 
between the cluster SZ peaks and the identified cluster cen- 
ters that has a uniform random distribution between 0 and 1.5 
Mpc, analogous to Figure 6. We also convolve both maps with 
ACT beams and a 5 ' Gaussian beam to approximate Planck 
resolution. The results of the extracted Y 5 qq values are shown 
as green diamonds in Figure 9. We also show the case where 
the simulated clusters are embedded in simulated CMB plus 
white noise maps, using white noise levels that are similar 
to Planck noise levels. These results are shown as purple 
hexagons in Figure 9. 

Planck- like noise is nearly white at these frequencies, so 
there is not much difference between the Planck-Wkc case and 

30 We doubled the number of maps and thus the cluster sample to shrink 
the error bars. This results in 984 embedded clusters in total. 



FIG. 10. — Shown is the filter response for differing beams and noise mod- 
els. The blue solid curve is a 5 ; Gaussian profile, which roughly approximates 
the filter response for a cluster profile convolved with the Planck beam since 
Planck's noise is close to white (see Eq. 1). The red dashed curve shows a 
1.4' Gaussian profile, which has been Fourier transformed and has had all 
Fourier modes set to zero for l < 2000, roughly analogous to the effect of the 
ACT noise. The profile was then Fourier transformed back. The green dotted 
curve shows the same as the red curve except for a 5' Gaussian profile. 


31 The instrumental noise in the ACT 218 GHz map is not low enough to 
remove all the primary microwave background signal from the 148 GHz map. 

32 Note that estimating the noise power spectrum as CMB plus white noise 
when filtering the ACT data results in very large errorbars since the matched 
filter is no longer optimal. High-pass filtering die ACT maps enough to mini- 
mize the effect of the red noise causes ringing of the true cluster signal, which 
also results in a low measurement if the signal is extracted from a position 
offset from the center. 
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It is interesting to see from Figure 9 that the amount of 
offset modeled in Figure 6 to match the ACT data, can also 
explain the discrepancy found by Planck . If this amount of 
miscentering actually exists, then it would be the sole expla- 
nation of the discrepancies. However, this amount of offset 
between BCGs and SZ peaks is much larger than the distri- 
bution shown in Figure 2. So it may be that the SZ signal 
is intrinsically low by some amount or that the optical weak- 
lensing mass calibration is biased high. It is also possible that 
the fraction of false detections in the optically-selected sam- 
ple is larger than 10%, or some non-zero amount of radio or 
infrared galaxy contamination is filling in the SZ decrements 
(as discussed in Section 5.4). If so, these effects would serve 
to lower the measured SZ signal by the same amount for both 
ACT and Planck. In that case, a smaller amount of miscen- 
tering would be required to match the remaining discrepancy 
found by these instruments. It is thus likely that a combina- 
tion of effects may be at work resulting in the measurements 
found by ACT and Planck. 

7. DISCUSSION 

From the analysis of simulated clusters embedded in ACT 
data presented in Figure 1, we expect robust recovery of the 
SZ flux from optically-selected MaxBCG clusters. This is as- 
suming that the position of the center of the gas concentra- 
tion is known. From the analysis of 52 MCXC clusters that 
fall within the ACT Equatorial and Southern survey regions 
shown in Figure 4, we find agreement between the expected 
and measured Y 5 oo c -Msqoc relation. Such agreement is con- 
sistent with that found in Planck Collaboration et al. (201 lb). 
Both of these figures taken together give confidence in the SZ 
flux recovery pipeline and the analysis of the ACT maps. 

In Figure 5, we find that the recovered ACT SZ flux from 
474 optically-selected MaxBCG clusters is lower than both 
the model expectations and the measured Planck values given 
in Planck Collaboration et al. (2011c). Since we expect from 
Figure 2 some offset between the positions of the BCGs and 
the centers of the gas concentrations, we explore this possibil- 
ity in more detail. Given the difference in resolution and noise 
properties between ACT and Planck , such an offset would re- 
sult in different measured SZ signals between the two instru- 
ments (see Figure 10). 

Figure 2 gives the offset distribution for clusters in common 
between the full MaxBCG and MCXC cluster catalogs. Using 
simulations, Figure 3 shows that this amount of offset is not 
enough to explain the measured SZ signal shown in Figure 5. 
Modeling the offset using the Johnston et al. (2007b) distribu- 
tion, which only includes misidentification of the BCG by the 
optical cluster finding algorithm, produces a slightly smaller 
amount of discrepancy between expected and measured SZ 
signals than what is shown in Figure 3. To match the mea- 
sured SZ signal found by ACT requires an offset distribution 
that is significantly larger. A distribution that is uniformly ran- 
dom between 0 and 1 .5 Mpc gives a better fit to the ACT data, 
as shown in Figure 6. This amount of offset also can explain 
the Planck measured discrepancy as shown in Figure 9. 

It is possible that the subsample of clusters in common be- 
tween MCXC and MaxBCG catalogs is not representative 
of the full MaxBCG cluster sample. This subsample is in 
fact special in that these clusters are found using both optical 
and X-rav selection techniques. When Planck measured the 
SZ signal for such a subsample, no significant discrepancy 
between the measured and expected SZ signals was found 
(Planck Collaboration et al. 201 lc). This is in contrast to their 


results for the full optically-selected sample. This suggests 
that the subsample of clusters in common to both catalogs is 
not representative of the full optically-selected sample, at least 
in some regard. 

We also assume throughout this work that the Amaud et al. 
(2010) profile accurately describes the gas profiles of clusters 
in the MCXC and MaxBCG catalogs. The good agreement we 
find in Figure 4 and also shown in Planck Collaboration et al. 
(2011b) suggests that this profile works well for clusters in 
the MCXC sample. However, this profile may not hold for 
the optically-selected MaxBCG sample. In particular, to mea- 
sure a lower SZ signal with ACT than found by Planck would 
require that the gas distribution is wider (more spread out) 
than that given by the Arnaud et al. (2010) profile. In such a 
case, given the finer resolution of ACT, it is conceivable that 
ACT could measure a lower signal than Planck when the SZ 
signal is extracted using an incorrect narrower profile for both 
datasets. However, since we expect clusters to be in hydro- 
static equilibrium to first order, the consequence of a wider 
gas distribution is a lower total mass within Msoq c if the gas 
mass fraction is fixed. This would imply lower normalizations 
of the L x ~M 5 qqc and A^oom “Mooc relations. That in turn 
would result in a lower a% value obtained from the MaxBCG 
sample than reported in Rozo et al. (2010), the latter of which 
is currently in agreement with <r 8 values obtained from X-ray- 
selected cluster samples (Henry et al. 2009; Vikhlinin et al. 
2009; Mantz et al. 2010). 

Some other recent papers have added more data and anal- 
ysis shedding light on the measured discrepancy in SZ signal 
found by Planck. Biesiadzinski et al. (2012) argue that there 
is not really a discrepancy if one considers the 2 a uncertainty 
on the N 2 oQ m -M 5 ooc relation in Rozo et al. (2009). However, 
if Rozo et al. (2010) used an N 2 oom ~ Mooc relation that had a 
normalization that is 2a away from the true value, then again 
it is surprising that they find a constraint on a* with an op- 
tical cluster sample that is in such good agreement with the 
a% constraints found using X-ray-selected samples. Such a 
shift in the N 2 oom - M 5 ooc normalization would also put the 
derived cosmological constraints in some tension with previ- 
ous results (Rozoetal. 2012a). Angulo etal. (2012), using 
Millennium-XXL simulations, argue along the same lines as 
Biesiadzinski et al. (2012) suggesting that the normalization 
of the richness-mass relation is incorrect (see last paragraph 
of Section 4 in that work). 

Bauer et al. (2012) tested the A^oom-A/sooc relation by veri- 
fying the weak-lensing cluster mass estimates upon which this 
relation was calibrated. They did this by using gravitational 
lensing magnification of type I quasars in the backgrounds of 
these clusters. Their results support the mass normalization 
of the Rozo et al. (2009) A^oom - A/rooe relation. Recent work 
by Rozo et al. (2012d) suggests that part of the Planck mea- 
sured discrepancy may be due to different X-ray instrument 
calibrations between Chandra and XMM-Newton and/or dif- 
fering analysis procedures when dealing with data from each 
instrument (see Figure 3 of that work). Rozo et al. (201 2d) 
argue that X-ray scaling relations derived from Chandra ob- 
servations may be systematically high compared to those from 
XMM-Newton observations. Since the A^oom - A/ 500 c relation 
was calibrated using an L x - A/ 5 0 Cc prior derived from Chan- 
dra observations (Vikhlinin et al. 2009), and the Amaud et al. 
(2010) profile was derived using XMM-Newton observations, 
this may be responsible for part of the discrepancy. Rozo et al. 
(2012a,b,c) include further discussion in this direction and 
claim that three ingredients may be responsible for the ten- 
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sion between measured and expected SZ signals measured 
by Planck for the MaxBCG sample. These ingredients are a 
mass-richness relation that is normalized high by 1 0%, Chan- 
dra derived X-ray masses from Vikhlinin et al. (2009) that are 
biased low by 20% due to non-thermal pressure, and XMM- 
Newton X-ray derived masses that are ~ 10 to 20% lower than 
those derived from Chandra. More data will clarify this issue. 

Further investigation is required using multi-wavelength 
techniques to fully contextualize the SZ signal measured by 
both Planck and ACT for this sample of optically-selected 
clusters. High-resolution SZ instruments such as CARMA 
and MUSTANG2 as well as more sensitive datasets such as 
those forthcoming from ACTPol and SPTPol will provide 
more clarity. This will yield a better understanding of the 
astrophysics of galaxy clusters. Such knowledge will also be 
vital to assess the use of galaxy clusters as probes of structure 
growth and the cosmological parameters of the Universe. 
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